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The DOUBLE-PRECISION program lets you obtain 15-digit 
accuracy with sine, cosine, arctangent, natural logarithm, 
exponential, and square root functions. These functions, 
together with LEVEL II BASICS "DEFDBL", and "DEFINT", provide 
a high degree of precision, for any applications that need a high 
level of accuracy. 

The compact subroutines are very efficient. If you own a LEVEL 
II - 4K machine, you can use these routines without concern. 
Enough memory will normally be left for your main program. 



CUSTOM MANUFACTURED IN U.S.A. BY RADIO SHACKS3 A DIVISION OF TANDY CORPORATION 



INTRODUCTION 

The DOUBLE -PRECISION program lets you obtain 15-digit 
accuracy with sine, cosine, arctangent, natural logarithm, ex- 
ponential, and square root functions. These functions, together 
with LEVEL II BASIC'S "DEFDBL", and "DEFINT", provide a 
high degree of precision, for any applications that need a high 
level of accuracy. 

The compact subroutines are very efficient. If you own a 
LEVEL II - 4K machine, you can use these routines without con- 
cern. Enough memory will normally be left for your main 
program. 

The three programs included are; 
SUBROUTINES 

The actual Double-Precision Subroutines are designed to be 
included in your program. To use them, simply execute a 
GOSUB to the particular routine required. Double-Precision 
Derived Functions may be included with slight modifications. 
You can DELETE any routines you won't need to conserve 
memory. 

ACCURACY TEST 1.1 

The Accuracy Test is a program designed to demonstrate the 
degree of precision that can be obtained with the Subroutine 
Functions. Each Accuracy test asks for a value to test, and an 
increment value. Each of the tests is then run 13 times, using 
your values. The results are graphically displayed, to aid you in 
determining the ACTUAL degree of precision obtained for any 
given number. In this manner, you can "Pre-test" some of the 
values that you plan to use in a long complex chain calculation 
to predict the final degree of accuracy, and the probable amount 
of error likely to occur. 

DOUBLE-PRECISION CALCULATOR 1.1 

Converts the TRS-80 to a calculator mode to give you the 
opportunity to easily experiment with the Double-Precision 
functions in degrees/angles. A unique "XXILLEGAL ARGU- 
MENT XX" sign will appear if the argument you enter will 
cause an error condition. 

All of the subroutines use the "Taylor Expansion Method" to 
support their precision. Extensive testing indicates that these 
routines will be accurate to 15 and 16 places in the majority. In 
a few tests we have determined that a 14th place digit could be in 
error by one or two units under certain conditions. 
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THEORY AND OPERATION 

SUBROUTINES 

All angles are presented in radians. 

The square root function is calculated by using the single- 
precision SQR function. The result is refined into double- 
precision, by dividing it into the argument, and then averaging 
with the quotient. This procedure is then repeated. With all other 
functions, the argument is reduced repeatedly, until it can be 
evaluated by a power series. The result is obtained by appro- 
priate expansion of the value (of the power series) the same 
number of times the argument was reduced. 



DESCRIPTION 



CALL 


FUNCTION 


ARGUMENT 


RESULT 


GOSUB 40100 


Natural Logarithm 


Z 


Z2 


GOSUB 40200 


Exponential 


Z 


Z2 


GOSUB 40300 


Cosine 


Z 


Z2 


GOSUB 40400 


Sine 


Z 


Z2 


GOSUB 40500 


Arctangent 


Z 


Z2 


GOSUB 40600 


Square Root 


z 


Z3 


GOSUB 40620 


PI/2 


none 


Z9 


GOSUB 40640 


Power Expansion 


(subroutine for 








above functions) 



IMPORTANT NOTE: 

In order to use these subroutines, you must define "Z" as a 
double-precision "variable" early in the program. (See "Level 
II" Manual, Page 4/3) At the start of your program, insert 
"DEFDBL Z". You must also keep in mind the fact that these sub- 
routines have some Variable Names reserved "For Subroutine 
Use Only". The Reserved Variable Names are; Z2, Z3, Z4, Z9, 1, 
10, II, 12, and 13. YOU CANNOT USE THESE NAMES IN YOUR 
PROGRAM, AND STILL MAINTAIN 16 DIGIT ACCURACY. 
Once entered, the argument in Z is destroyed. Most functions 
output by the variable Z2. SQR outputs by the variable Z3. PI/2 
outputs as Z9. You must not redefine the value held in Z9. Once 
calculated, the subroutines for the trigonometric functions ex- 
pect to find the value for PI/ 2 stored in Z9. 



SUBROUTINE 
SPECIFICATIONS 

NATURAL LOGARITHM: 

Purpose of Function Computes Z2=L0G(Z) 

Input Z 

Output Z2 

Reserved Variables Used 1, 10, 12, 13 

Error Conditions Z< =0 

Other Subroutines Used 40600, 40640 

Calling Procedure GOSUB 40400 

EXPONENTIAL: 

Purpose of Function Computes Z2=EXP(Z) 

Input Z 

Output Z2 

Reserved Variables Used I, 10 

Error Conditions overflow for Z>88 or Z488 

Other Subroutines Used none 

Calling Procedure GOSUB 40200 

COSINE: 

Purpose of Function Computes Z2=COS(Z) 

Input Z (in Radians) 

Output Z2 (in Radians) 

Reserved Variables Used Z9 

Error Conditions none 

Other Subroutines Used 40620, 40400 

Calling Procedure GOSUB 40300 

SINE: 

Purpose of Function Computes Z2=SIN(Z) 

Input Z (in Radians) 

Output Z2 (in Radians) 

Reserved Variables Used I, 10, 12, Z9 

Error Conditions none 

Other Subroutines Used 40620 

Calling Procedure GOSUB 40400 

ARCTANGENT: 

Purpose of Function Computes Z2=ATN (Z) 

Input Z (in Radians) 

Output Z2 (in Radians) 

Reserved Variables Used I, 10, II, 12, 13, Z3, Z9 

Error Conditions none 

Other Subroutines Used 40620, 40640 

Calling Procedure GOSUB 40500 



SQUARE ROOT: 

Purpose of Function Computes Z3=SQR(Z) 

Input Z 

Output Z3 

Reserved Variables Used none 

Error Conditions Z<=0 

Other Subroutines Used none 

Calling Procedure GOSUB 40600 

PI/2: 

Purpose of Function Computes TT/2 

Input none 

Output Z9 

Reserved Variables Used Z2, Z4 

Error Conditions none 

Other Subroutines Used 40500 

Calling Procedure GOSUB 40620 



ACCURACY TEST 

This program performs four tests to demonstrate the accuracy 
of the double-precision functions. The results show the maxi- 
mum degree of error allowed for any given number: 



TEST 

SIN»+COS«-l 
ATN(SIN/COS) 
LOG(EXP) 
(SQR) 2 



BEST POSSIBLE RESULT 


Results = Argument 
Results = Arguments 
Results = Arguments 



CALCULATOR 

The Calculator program lets you enter any function of the 
form; 

XXX(Argument) 

wherein XXX is SIN, COS, ATN, LOG, EXP, or SQR. The argu- 
ment cannot be an expression. All angles are in degrees. If the 
argument causes an error condition, an "XXILLEGAL ARGU- 
MENT XX" message appears. If XXX is not one of the six defined 
functions, a "WHAT?_" will appear. Type END to stop the 
program. 



EXAMPLES; 






Computer; 


You type; 


Computer displays; 


? 


SIN(30) 


.4999999999999999 


?_ 


COS(30) 


.8660254037844387 


? 


LOG(l) 





? 


EXP(l) 


2.718281828459043 


? 


ATN(l) 


45 


? 


SQR(-l) 


XX ILLEGAL ARGUMENTXX, 


? 


END 


READY > _ 



DERIVED FUNCTIONS 

You can also achieve Double-Precision results for Derived 
Functions — by simply calling out "DEFDBL" for all Constants, 
and the appropriate GOSUB for each part of the Derived 
Function that requires the Double Precision Subroutine Series. 
We have listed the Derived Functions, and the function as 
expressed in terms of LEVEL II Basic; 



FUNCTION 



FUNCTION EXPRESSED IN 
TERMS OF LEVEL II BASIC 
FUNCTIONS 



SECANT 
COSECANT 
COTANGENT 
INVERSE SINE 

INVERSE COSINE 

INVERSE SECANT 

INVERSE COSECANT 

INVERSE COTANGENT 
HYPERBOLIC SINE 
HYPOBOLIC COSINE 
HYPERBOLIC TANGENT 

HYPERBOLIC SECANT 
HYPERBOLIC 

COSECANT 
HYPERBOLIC 

COTANGENT 

INVERSE HYPERBOLIC 
SINE 

INVERSE HYPERBOLIC 
COSINE 

INVERSE HYPERBOLIC 

TANGENT 

INVERSE HYPERBOLIC 
SECANT 

INVERSE HYPERBOLIC 
COSECANT 

INVERSE HYPERBOLIC 
COTANGENT 



SEC(X) = l/COS(X) 
CSC(X) = 1/SIN(X) 
COT(X) = 1/TAN(X) 
ARCSIN(X) = ATN(X/SQR 

(— X*X+1)) 
ARCCOS(X) = -ATN(X/SQR 

(— X*X+1))+1.5708 
ARCSEC(X) = ATN(SQR(X*X— 1)) 

+(SGN(X)— 1)*1.5708 
ARCCSC(X) = ATN(1/SQR 

(X*X-1))+(SGN(X)-1)*1.5708 
ARCCOT(X) = — ATN(X)+ 1.5708 
SINH(X) = (EXP(X)-EXP(-X))/2 
COSH(X) = (EXP(X)/EXP(-X))/2 
TANH(X) = -EXP(-X)/(EXP(X) 

+EXP(— X)*2+l 
SECH(X) = 2/(EXP(X)+EXP(— X)) 

CSCH(X) = 2/(EXP(X)— EXP(— X)) 

COTH(X) = EXP(— X)/(EXP(X) 
— EXP(— X))*2/l 

ARGSINH(X) = LOG(X+SQR 
(X*X+1)) 

ARGCOSH(X) = LOG(X+SQR 
(X*X-1)) 

ARGTANH(X) = LOG ((1+X)/ 
(l-X))/2 

ARGSECH(X) = LOG((SQR(— X*X+ 
1)+1)/X) 

ARGSCH(X) = LOG((SGN(X)*SQR 
(X*X+1)+1)/X) 

ARGCOTH(X) = LOG((X+l)/ 
(X-l))/2 
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DOUBLE PRECISION F'RGE 



40606 END 

40020 - DOUBLE PRECISION ROUTINES 

40040 -'COPYRIGHT TANDY CORP. 1978 

40050 '• 

40052 ■' *** NATURAL LOG *** 

40054 ' 

40100 22=LOG<Z>:I0=0:I2=SGN<Z2):IF I2<0 Z=l/Z 

40110 IF Z<1. 065 THEN 40130 

4012O GOSUB 40600 :Z=Z3: 10=10+1: GOTO 40110 

40130 Z=<Z-1>/(Z+1) :Z2=Z*Z:Z3=Z 

40140 I 3=9: GOSUB 40640 :Z=Z+Z 

40150 IF 10=0 THEN 48170 

4016O FOR 1=1 TO IO:Z=Z+Z:NEXT 

40170 Z2=Z*I2: RETURN 

40180 ' 

40182 ■' *** EXPONENTIAL *** 

40184 ' 

40200 Z2=EXPCRBS<Z>) : I0=O 

40210 IF Z*ZC 004 THEN 40230 

40220 Z=Z/2: 10=10+1: GOTO 48210 

40230 Z2=1:F0R 12=8 TO 1 STEP -1:Z2=Z2*Z/I2+1:NEXT 

40240 IF 10=0 THEN 40268 

40250 FOR 1=1 TO I0:Z2=Z2*Z2:NEXT 

40260 RETURN 

40270 ■' 

40272 ' *+* COSINE *** 

40274 ' 

40300 GOSUB 40620 "COS 

40310 Z=Z9-Z 

40320 / 

40322 ' **+ SINE +** 

40324 '' 

40400 GOSUB 40620 'SIN 

48410 10=0: I2=SGN<Z) :Z2=Z9+Z9:Z=ftBS(Z) :Z=Z-INT<ZA.Z2+Z2>>*<Z2+Z2> 

40420 IF Z>Z2 THEN Z=Z-Z2: I2=-I2 

48430 IF Z>Z9 THEN Z=Z2-Z 

40440 IF flBSCZK 063 THEN 40460 

40450 Z=Z/3: 10=10+1: GOTO 48440 

40460 Z2=-Z*Z : Z2=< « CZ2/42+l>*Z2/20+l>*Z2/6+l>*Z 

40470 IF 10=0 THEN 40490 

40480 FOR 1=1 TO I0:Z2=(3-4*Z2*Z2)*Z2:NEXT 

40490 Z2=Z2*I2: RETURN 

40492 " 

40494 '' *** ARCTANGENT *** 

40496 ■' 

40580 GOSUB 40628 

40518 1 0=0 : 1 1=0 : 1 2=SGN < Z ) 

40520 Z=ABS<Z>:IF Z>1 Z=1/Z:I1=1 

40530 IF Z<. 077 THEN 40560 

40540 Z2=Z:Z=Z*Z+1: GOSUB 48680 :Z=Z2A.Z3+1) 

40550 I 0=1 0+1 .GOTO 40530 

40560 Z3=Z:I3=-11:G0SUB 40640 
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DOUBLE PRECISION PFlGE 



49570 IF 10=0 GOTO 40578 

40574 FOR 1=1 TO I0:Z=Z+Z:HEXT 

48573 Z2=Z 

40582 IF 11=1 Z2=Z9-Z2 

40586 Z2=Z2*I2 

40590 RETURN 

40592 ' 

40594 '' *** SQUARE ROOT *** 

40596 " 

40600 Z3=SQR<Z> : Z3=<Z3+Z/Z3>/2 : Z3=(Zi+Z/Z3)/2 : RETURN 

46610 ' 

40612 ' *** Z9 = PI/2 *** 

40614 ' 

40620 IF Z9O0 RETURN 

40630 Z4=Z : Z=l : GOSUB 40510 : Z9=Z2+Z2 : Z=Z4 : RETURN 

40632 -' 

40634 ' *#* POWER EXPANSION SUBROUTINE *** 

40636 ' 

40646 Z2=Z*Z:F0R 1=3 TO ABSU3) STEP 2:Zi=SGNCI3>*Z3*Z2 

40650 Z=Z+Z3/ I : NEXT : RETURN 
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RCCURfiCV TEST PRGE 



100 ' DOUBLE PRECISION LEVEL II BASIC FUNCTIONS 

110 ' SIN, COS fiTN, LOG, EXP, 5QR 

120 -' 

130 ' ACCURACY TEST 

140 -COPYRIGHT TANDY CORP. 1978 

145 DEFDBL Z:DEFINT I-N:ONERRORGOTO9000 

150 CLS : PR I NTTAB ( 20 ) " ACCURACY TEST 1.1": PR I NT 

160 FGRJ=1T04:PRINTTAB(15);,T; "> "j : GOSUB300 : NEXT 

165 PR I NTTAB (15); " ) STOP": PRINT 

170 PR I NTTAB ( 15) ; : INPUT "ENTER SELECTION"; J 

ISO IFJ=0THEN ONERRORGOTOO : END ELSE GOSUB 20O:GOTO150 

280 H=0 : GOSUB40O : CLS : PR INTT AB ( 15 ) i : GOSUB30O : PR I NT " ARGUMENT " T AB ( 32 ) " RESULT 

210 ER=0 : ZA=ZS+CDBL ( M )*2I : Z=ZA 

220 ONJGOTO10O0J 2000.. 3800, 4000 

230 IFER=0THENPRINTZA; TAB(25); Z2 

240 M=M+1:IFM<13THEN 210 

250 PRINTG970, "PRESS ENTER TO CONTINUE".; : I NPUTA*: RETURN 

300 ONJGOTO310, 320.. 330, 346 

310 PRINT" (SIN t 2) + (COSt 2) -1":RETURN 

320 PRINT"ATN (SIN / COS)" .RETURN 

330 PR I NT "LOG (EXP)" .RETURN 

340 PR I NT " ( SQR ) t 2 " : RETURN 

4S0 PRINT: PR I NT "ENTER INITIAL ARGUMENT VALUE "TAB (32); .INPUT ZS 

410 PR I NT "ENTER INCREMENTING VALUE "TAB (32); : INPUTZI : RETURN 

1000 GOSUB40400:IFER=1THEH RESUME 230 

1010 ZB=Z2:Z=ZA:GOSUB40306:1FER=1THEN RESUME 230 

1020 Z2=ZB+ZB+Z2+Z2-1D0 : G0T0230 

2000 GOSUB4O620 : ZC= I NT ( ZA/Z9 ) : ZD=ZA-Z9*ZC : Z=ZD : GOSUB40400 : I FER=1THEN RESUME 230 

2010 ZB=Z2:Z=ZD:GOSUB40300:IFER=1THEN RESUME 230 

2020 Z=ZB/Z2:GOSUB40500:IFER=1THEN RESUME 230 

2030 Z2=Z2+Z9*ZC : G0T0230 

3000 GOSUB46200:IFER=1THEN RESUME 230 ELSE Z=Z2 

3010 GOSUB40100:IFEF.:=1THEN RESUME 230 ELSE GOTO 230 

4000 GOSUB40600:IFER=1THEN RESUME230 ELSE Z2=Z3*Z3 : G0T0230 

9000 PRINTZA;TAB(25)"+* ILLEGAL ARGUMENT **" :ER=1: RETURN 

40000 END 

40020 -' DOUBLE PRECISION ROUTINES 

40050 " 

40052 ■" *** NATURAL LOG *** 

40054 •" 

40100 Z2=LOG(Z):I0=0:I2=SGN(Z2):IF I2<0 Z=l/Z 

40110 IF Z<1. 065 THEN 40130 

40120 GOSUB 40600 :Z=Z3: 10=10+1: GOTO 4O110 

40130 Z=(Z-1)/(Z+1) :Z2=Z*Z:Z3=Z 

40140 I 3=9: GOSUB 40640 :Z=Z+Z 

40150 IF 10=0 THEN 40170 

40160 FOR 1=1 TO IO:Z=Z+Z:NEXT 

40170 Z2=Z+I2: RETURN 

40180 - 

40182 - **+ EXPONENTIAL *** 

40184 -' 

40200 Z2=EXP ( ABS ( Z ) ) : I 0=0 

40210 IF Z+Z<. 004 THEN 40230 

40220 Z=Z/2: 10=10+1: GOTO 40210 
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RCCURflCV TEST PAGE 



46230 Z2=1:F0R 12=8 TO 1 STEP -1:Z2=Z2*Z/I2+1:NEXT 

40240 IF 10=0 THEN 40260 

40250 FOR 1=1 TO I0:Z2=Z2*Z2:NEXT 

40260 RETURN 

40270 ' 

40272 '' *** COSINE *** 

40274 '' 

40300 GOSUB 40620 'COS 

40310 Z=Z9-Z 

40320 •' 

40322 " **+ SINE *** 

40324 -' 

40400 GOSUB 40620 "SIN 

40410 10=0: I2=SGN<:.Z) :Z2=Z9+Z9:Z=ABS<Z) :Z=Z-INT<ZAZ2+Z2))>KZ2+Z2) 

40420 IF Z>Z2 THEN Z=Z-Z2: I2=-I2 

40430 IF Z>Z9 THEN Z=Z2-Z 

40440 IF ABS<ZX. 063 THEN 40460 

40450 Z=Z/3: 10=10+1: GOTO 40440 

40460 Z2=-Z*Z : Z2=< (. <:Z2/42+l>+Z2/20+l>*Z2/6+l>*Z 

40470 IF 10=0 THEN 40490 

40480 FOR 1=1 TO I0:Z2=<3-4*Z2*Z2>*Z2:NEXT 

40490 Z2=Z2*I2:RETURN 

40492 ' 

40434 ' *** ARCTANGENT *** 

40496 ' 

40500 GOSUB 40620 

40510 1 0=0 : 1 1=0 : 1 2=SGN (. Z ) 

40520 Z=ABS<Z> : IF Z>1 Z=1/Z:I1=1 

40530 IF Z<. 077 THEN 40560 

40540 Z2=Z:Z=Z*Z+1: GOSUB 40600 :Z=Z2AZ3+1> 

40550 I 0=1 0+1 .GOTO 40530 

40560 Z3=Z:I3=-11:G0SUB 40640 

40570 IF 10=0 GOTO 40578 

40574 FOR 1=1 TO I0:Z=Z+Z:NEXT 

40578 Z2=Z 

40582 IF 11=1 Z2=Z9-Z2 

40586 Z2=Z2*I2 

40590 RETURN 

40592 ' 

40594 ' *** SQUARE ROOT ***■ 

40596 •' 

40600 Z3=SQR<Z> : Z3=<Z3+Z/Z3)/2 : Z3=<Z3+Z/Z3)/2 : RETURN 

40610 ' 

40612 ' *** Z9 = PI/2 *** 

48614 ' 

40620 IF Z9O0 RETURN 

40630 Z4=Z:Z=1: GOSUB 40510 : Z9=Z2+Z2 : Z=Z4 : RETURN 

40632 ' 

40634 •' *** POWER EXPANSION SUBROUTINE *** 

48636 ' 

40640 Z2=Z*Z:F0R 1=3 TO ABS<I3) STEP 2:Z3=SGN<I3)*Z3*Z2 

40650 Z=Z+Z3/I : NEXT : RETURN 
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CRLCULRTOR PROGRAM PAGE 



106 - DOUBLE PRECISION LEVEL II BASIC FUNCTIONS 

110 ■' SIN/ COS, ATN, LOG, EXP, SQR 

126 - 

ISO "• CALCULATOR 

140 - COPYRIGHT TANDY CORP. 1978 

150 DEFDBL Z.DEFINT I-N.ONERRORGOTO400 

160 F0RM=1T06:READ F$<M) : NEXT : GOSUB40628 

170 DRTR SIN, COS, HTN, LOG, EXP, SQR 

ISO CLS : PR I NTTRB < 15 ) " DOUBLE PREC I S I ON CRLCULRTOR 1. 1 " : PR I NT 

185 PRINTTAEKlyV'RNGLES RRE IN DEGREES" 

140 PRINT"ENTER FUNCTIONS — IE. 'SIN<M>" ~ 'END- TO STOP 

200 INPUT R*:F1*=LEFT$<A$, 3) : IFR$="END"ORR$="STOP"THENONERRORGOTO0:END 

210 F0RM=1T06 : IFF$ai>=Fl$THEN220ELSENEXT : PR I NTTRB < 15) "WHRT?" : GQTO200 

220 Z=CDBL':VRL(MID*<R*, 5, LEN<fi$>-5>>> 

220 IFH=1GRM=2THEN Z=Z+Z9/yDl 

240 ONMGOSUB40400, 40300, 40500, 48188, 48288, 406O0: IFM=6THENZ2=Z3 

250 I FM=3THENZ2=Z2+yDl/Zy 

260 PRINTTRB(15>; Z2:GOTO200 

400 PR I NTTRB CIS)"** ILLEGAL ARGUMENT **" :RESUHE20O 

40000 END 

40020 •■ DOUBLE PRECISION ROUTINES 

40050 ■'' 

40052 •'' *** NATURAL LOG *** 

40054 •'' 

40100 Z2=L0G<Z>: I0=0:I2=SGN<Z2>:IF I2<0 Z=l/Z 

40110 IF Z<1. 065 THEN 40130 

40120 GOSUB 48688 :Z=Z3: I 6= I 8+1: GOTO 40110 

40130 Z=<Z-1>/':Z+1>:Z2=Z+Z:Z3=Z 

40140 1 3=9: GOSUB 4O640:Z=Z+Z 

4015W IF 10=0 THEN 49178 

40160 FOR 1=1 TO IO:Z=Z+Z:NEXT 

4O170 Z2=Z*I2:RETURN 

40180 "' 

40182 " *** EXPONENTIAL *** 

40184 •' 

4O200 Z2=EXP(ABS(Z>>:I0=0 

40210 IF Z*Z<. 004 THEN 40230 

40220 Z=Z/2 : I 0= I 8+1 : GOTO 48218 

40230 22=1: FOR 12=8 TO 1 STEP -1:Z2=Z2+Z/I2+1:NEXT 

40240 IF 10=0 THEN 40260 

40250 FOR 1=1 TO I8:Z2=Z2*Z2:NEXT 

40260 RETURN 

40270 " 

40272 -' *** COSINE *** 

40274 -' 

40300 GOSUB 40620 "COS 

40310 Z=Z9-Z 

40320 -' 

48322 ' *** SINE *** 

40324 ■' 
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CRLCULflTOR PROGRAM PAGE 



46406 GOSUB 40620 "SIN 

40410 I0=0:I2=SGN<Z>:Z2=Z3+Z9:Z=fiBS<Z>:Z=Z-INT<ZA22+Z2>>*<:Z2+Z2> 

40420 IF Z>Z2 THEN Z=Z-Z2: I2--I2 

40430 IF Z>Z9 THEN Z=Z2-Z 

40440 IF RBS<2K. 063 THEN 40460 

40450 Z=Z/3: 10=10+1: GOTO 40440 

40460 Z2=-Z*Z : Z2=< < <Z2/42+l>*Z2/20+l>*Z2/6+l>*Z 

40470 IF 10=0 THEN 40490 

40480 FOR 1=1 TO I0:Z2=<3-4*Z2*Z2>*Z2:NEXT 

40490 Z2=Z2*I2: RETURN 

40492 •' 

40494 ' *++ ARCTANGENT *** 

40496 ' 

40500 GOSUB 40620 

40510 1 0=0 : 1 1=0 : 1 2=SGN ( Z > 

40520 Z=ABS<Z>:IF Z>1 Z=1/Z:I1=1 

40530 IF Z<. 077 THEN 40560 

40540 Z2=Z : Z=Z*Z+1 : GOSUB 40600 : Z=Z2 A Z3+1 ) 

40550 I 0=1 0+1 .GOTO 40530 

40560 Z3=Z:I3=-11: GOSUB 40640 

40570 IF 10=6 GOTO 40578 

40574 FOR 1=1 TO IO:Z=Z+Z:NEXT 

40578 Z2=Z 

40582 IF 11=1 Z2=Z9-Z2 

40586 Z2=Z2*I2 

40590 RETURN 

40592 * 

40594 ' *** SQUARE ROOT *** 

40596 ' 

40600 Z3=SQR C2> : Z3= < Z3+Z/Z3 > /2 : Z3= < Z3+Z/Z3 > /2 : RETURN 

40610 "' 

40612 ' *** Z9 = PI/2 *** 

40614 ' 

40620 IF Z9O0 RETURN 

40630 Z4=Z:Z=1: GOSUB 40510 : Z9=Z2+Z2 : Z=Z4 : RETURN 

40632 ' 

48634 ' *** POWER EXPANSION SUBROUTINE *** 

40636 ' 

40640 Z2=Z*Z:F0R 1=3 TO ABSU3) STEP 2:Z3=SGN<I3)*Z3*Z2 

40650 Z=Z+Z3/ I: NEXT: RETURN 
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IMPORTANT NOTICE 

ALL RADIO SHACK COMPUTER PROGRAMS ARE DISTRIBUTED ON AN 
"AS IS" BASIS WITHOUT WARRANTY 

Radio Shack shall have no liability or responsibility to customer or any other 
person or entity with respect to any liability, loss or damage caused or alleged to 
be caused directly or indirectly by computer equipment or programs sold by 
Radio Shack, including but not limited to any interruption of service, loss of 
buisness or anticipatory profits or consequential damages resulting from the use 
or operation of such computer or computer programs. 

NOTE: Good data processing procedure dictates that the user test the program, 
run and test sample sets of data, and run the system in parallel with the 
system previously in use for a period of time adequate to insure that 
results of operation of the computer or program are satisfactory. 
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